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In the framework of the quasipotential method we calculate the contributions 
of the kind (Za) e m1/m^ to the energy spectrum of the muonium n 3 Si states. 
Numerical value of obtained correction for 2 3 S\ -j- l 3 Si muonium fine structure 
interval is equal to 0.19 MHz. 
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The investigation of muonium and positronium fine structure represents one of the 
basic tests of quantum electrodynamics, which is sensitive to radiative corrections of 
higher order on a [ffl. Many papers 0, |3|, [|, [|] are devoted to the calculation of different 
contributions to the fine structure of hydrogen-like atom energy levels. The interest 
to this problem remains unchanged [0, [7], ||. The progress, achieved in the last years 
during the calculation of logarithmic contributions of order a 6 In a in the positronium fine 



structure intervals (2 3 5i 4- l 3 Si, 2 3 S'i -4 2 3 Pj) |n| [U|, [I2|, doesn't abolish the necessity of 
calculation of higher order corrections 0(a 6 ) As discussed in the paper 0, there is 
some difference in the calculation of the corrections (Za) 6 mf/m 2 , obtained in |8], [I]]. The 
development of experimental methods, based on Doppler-free two-photon spectroscopy, 
allows the "large" structure intervals to be measured for the muonium and the positronium 



A7 ™ 9 3o . 1 3 a \ / 1233607218, 9 ± 10, 7 MHz 
AE Ps (iii.l b!) - j 1233607216) 4 ± 3 2 MRz (1) 

AE e ^(2 3 S 1 4 l 3 Si) = 2455527936 ± 120 ± 140 MHz. (2) 

The frequency of Doppler-free two-photon transition 15* -j- 25* in the hydrogen atom as 
well as hyperfine splitting of the ground state of hydrogen atom represents the quantity, 
which was measured with high accuracy JjlJ. The increase of the experimental accuracy 
of the muonium fine structure interval measurements (just as for positronium), planned 
in the near future, makes very actual the calculation of radiative corrections of higher 
order on a. In this paper we have performed studies of the recoil contributions of order 
0(a 6 ) in the muonium fine structure. The contribution of such order to the muonium 
hyperfine structure was obtained in |17| . There exist many approaches for description 
of relativistic energy spectrums of two-particle bound states in quantum electrodynamics 
H [|, [yl [19|, p0| . They all differ in organization of the calculation: the bound state 



equation for two-particle system, the construction of the interaction operator, the degree 
of complication in the calculation of the definite order corrections in the energy levels. 
But all such methods give equivalent results in the fixed order of perturbation theory on 
the small parameter a and mi/m 2 . Our calculations are based on the Schrodinger-type 
local quasipotential equation pi 



^ ^ i>M(p) = I 7^V(p, q, M)ij M (q), (3) 



\2fi R 2jj, R J J (2tt 

where b 2 = Ef — mf — E\ — m 2 ., fiR = ExE^/M is relativistic reduced mass, M = E\ + E2 
is the bound state mass, mi, 777-2 are the masses of the electron and the muon. As an 
initial approximation of quasipotential V(p, q, M) for the bound state system (e~fi + ) we 
choose the ordinary Coulomb potential. On the basis of equation (3) in |22[] we have 
obtained some relativistic corrections ma 6 in the positronium fine structure from the 
one-photon, two-photon interaction and the second order perturbation theory. First of 
all, the contribution of order [Zafmxjvn^ 

2777277° 



2 



appears from the condition of quantization of the energy levels for Coulomb interaction 

(5) 

which we have transformed for binding energy B. 



2 9 ' 

/4 n 2 



1 One-photon interaction contribution to the fine 
structure 

The basic contribution to the energy spectrum of two-particle bound state (/i + e~) is 
determined by one-photon interaction. The construction of one-photon quasipotential in 
the system of two spinor particles was done in J22]. . It is useful to note, that procedure 
of quasipotential derivation may be essentially simplified in this case, just as for two- 
photon and three-photon exchange diagrams by introducing the relativistic operator of 
projection onto 3 Si muonium state: 

^_ 1 (pi + mi)(l +7o)g(~P2 + m 2 ) 
2^2 \fe\ + m 1 ^e 2 + m 2 

(where pi, p 2 are four-momenta of electron and muon in the initial state, is the muonium 
polarization vector. Doing the decomposition of all relativistic factors over two small 
parameters: |p|/m 12 (\p\ ~ Za is the momentum of relative motion) and m\jm 2 in order 
to extract the contribution of six order over a and the first order over m\jm 2 we can 
present the interaction operator in coordinate space in the following manner: 



Vl( r) = - *i - (l+ - ^ { ?V) (l + - (7, 

lm\r A \ m 2 J 4mfr d V m 2 



-aZa ... Zab 2 . . . 

-5(r) = V c (r) + AVi(r) 



"im\m 2 m\m 2 r 

It is necessary to point out, that the part of quasipotential (7) AVi taking together 
with quantization condition (5) reproduce correctly the well-known energy spectrum of 
muonium S-wave states with the accuracy 0(a 4 ) |TB|, ^U|. The terms of AVi m our 
approach give also the corrections 0(a 6 ) in the energy levels. The reason is that the 
relativistic reduced mass {Ir and h 2 are dependant from a. Averaging (7) over Coulomb 



wave functions [24] and extracting necessary order contributions, we have obtained: 

3ml(Zaf 



4m 2 n 5 V n 

The quasipotential of one-photon interaction contains also a number of other terms, 
which lead to the energy corrections of order OiZaf '. They may be found when con- 
structing Vi 7 with the accuracy to terms of fourth order over \p\/mi, 2 , \q\/mi )2 in the 
form 

AV 2 (p, q,M = {—- 4 3 -) + 3 + — " 9 

k 2 L 16mf V m 2 ) 9omf \ m 2 J 
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if + t)(P0) f Q + 13mA _ (p 2 + cF)b 2 / _ mA _ (pq)b 2 / + lm x 



96mf \ m 2 / 96mf V va^J 48mf V m 2 

Let consider calculation of O(Za) 6 corrections from AV2. A series of AV2 terms will 
cause the divergent integrals in the energy spectrum. The reason behind this divergence 
is that, in deriving (9), we expanded all quantities determining the quasipotential in the 
relative momenta, which are proportional to a. A typical divergent integral with respect 
to momenta has the form / p 2 dpil) n s{p)- The relativistic correction of order a 6 in this case 
is determined by the residue of the integrand at the pole of the wave function ip n s(p) [P~8" 



The result of this integral calculation for arbitrary principal quantum number n has the 
form: 

dp f [3 + 2(n-l)(n + l)l 9 , 

Wfk^ 1 = — * — la M = 0) - (10) 

Taking into account (9) under averaging AV 2 we have calculated the relativistic corrections 
of necessary order for the energy levels with arbitrary principal quantum number n: 

m 2 {Zaf ( 139 171n2 73 1 43 1 17, iXn 1 r ^ ... \ 
where ip(z) = d\nT(z)/dz, C=0. 5772156649... is the Euler constant. 



2 Second order of perturbation theory 

In the second order of perturbation theory, the correction to the muonium energy spectrum 
is determined in our approach by the sum of two terms ||25|| : 



^ =< ,,|av,i,, >< ci^i* > + i ~ > 



dB 



k=l,k^n 



B c n ~Bl 



(12) 

The quasipotential (9) explicitly depends on the bound state energy B (the factors b 2 and 
Hr). Bearing in mind that to a precision adopted here, the relation db 2 /dB = 2/x holds, 
we obtain: 



A£ 4 =< CI A Wn >< 



' dB 



1 



m 2 



TV 



(13) 



The sum over Coulomb states in (12) is defined by the reduced nonrelativistic Green's 



function |L9], JZQ, |2^, |28|, |29], |30], |3l| , which has the following partial expansion: 

G n (r,r *,B) = £ g nl (r, r', B)Y lrn {n)Y* m {n') (14) 



Lm 



The radial function g n i(r,r', B n ) was obtained in |3T| as an expansion over Laguerre's 
polynomials. In the case of S-wave states we have for it the next expression: 



9no(r, r',B n ) 



4Zafi 2 



n 



x-\-x 

e ~ 



m=l,m^n 



m{m — n) 



(15) 
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where x = 2fiZar/n, L™ are the ordinary Laguerre's polynomials, which may be found 
by means the following relation: 

LZ(x) = ^ r [±) (e-^" +m ). (16) 

The reduced Coulomb Green's function (RCGF) depends from two variables r and r'. 
Due to the appearance of the 5 - like potential in (9) the RCGF is needed for r = 0. 
The necessary expression for RCGF was obtained by means of Hostler's representation 
for Coulomb Green's function in |32]| , after subtraction of pole-like term: 

^^..v.-.gfcf^ (17) 

nirx ~^ si [(n — s)\\ z 

{, xr / \ , , 2(n — s) + 3 — x , . ] 
Un-s) [ip(n + 1) - 2ip(n - s + 1) ^ h In x\ + 1 > , 



where ip(z) = d\nT(z)/dz. Contrary to the paper |32| this formula doesn't contain free 
two-particle Green's function G^(r) = —fiRe~ ZatlHr /27ir, which controls the iteration part 
of quasipotential. Its contribution to the energy spectrum will be find separately. Let 
consider, for example, calculation of energy correction in the second order of perturbation 
theory, depending from 5-like part of the quasipotential and the term of AV\ ~ 1/r 2 . 
This contribution can be represented in the form: 

5B _ ^(Zaf " fc n! 

3mim2n 4 ^ [n — k)\(k\) 2 



OCj 



{2k 4- 3 — r I 
k[i){n + 1) - 2^{k + 1) + lnx] + 1 . 

The expression (18) contains the integrals of two types with the power and the loga- 
rithmic functions correspondingly. Calculation of the first integral over x variable leads 
to the result 

h ~Jo X 6 L ^ [X) - ( n -l)\T(2-k) ' li9j 
So, we see that the sum over k contains only one term with k=l. Second integral of (18) 

J 2 = [°° x k -i \nxe- x L l n _ x {x) = (2 ~ k ^~^ ty( k ) + if) (2 — k) — ip(n + 1 - k)} (20) 
Jo [n — ly. 

gives rise the sum of the kind 

» fc V>(2-fc)(2-fc) w - 1 n-l 

ti } (n-k)\k\ " ~ + C ' (21) 
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where C = lim r , 



Inn + J2m=i ~] = 0.57721566.... Taking into account that 
r(n - k) 



lim 



= (-!)"->-!)!, 



fc-m r(i - k) 

we have obtained conclusively respective correction (18): 

TH 2 1 

5B = --±(Za) 6 -(4n + 9). 
bin? n 4 



(22) 



(23) 



Table 1. Second order contributions of the perturbation theory, defined by 
the RCGF and quasipotential (9) in the units ^ Za ^ mi 



AV 1 



Zab 2 



(Za) 2 fi R q _|_ 2mi 



Za(rV) , 
4mfr 3 



I —J— 4m l 

m2 



3mim2 



5(f) 



mim2r 



2n 2 -5n+l 
4n 6 



(^)Vii /i I 2mi \ 
2m 2 r 2 ^ ^ m 2 > 



2-n 

r)5 



2n+3 
2re 4 



n 2 +3w-l 
4n 5 



(4n+9) 
6n 4 



Za(rV) q _|_ 4m i ' 



2n 2 -5ra+l 
4n 6 



n 2 +3n-l 
4n 5 



n(n-l)(n+l) 
24n 6 



(n 2 -6n+8) 
12n 5 



3mim2 



5(f) 



(4n+9) 
6n 4 



(n 2 -6ra+8) 
12n 5 



Similar calculations may be carried out to find the contributions of the other terms of 
quasipotential (9) in the second order of perturbation theory. The respective results of 
such calculation represented in the Table 1. The full correction, connected with reduced 
Coulomb Green's function without considering G* in the energy spectrum is: 

25 _ 3 _49_ _3_\ m\(Zaf 
24 ~ n ~ 



A Ik 



(24) 

24n 2 2n 3 J m 2 n 3 V J 

Let consider now omitting in the RCGF contribution of free two-particle propagator to 
the correction AB^ . It is convenient to perform necessary calculations in the momentum 
representation. Taking in mind that 

(27r) 3 5(p - q) 



Gf(p,q,M) = 



b 2 



(25) 



and J- like term of quasipotential (9) already contains the muon mass in the denominator, 
we may represent the iteration-type correction as: 



AB, 



2fi(Za) 2 nij nS (0) 
3m 2 mf 



pk 



dq 



(27r) 3 (<f + l¥ 2 )' 



(26) 
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p 1 = (E 1 ,p) (E 1 + k ,k) q 1 = (E 1 ,q) 




p 2 = (E 2 , -p) (E 2 - k , -k) q 2 = (E 2 , -q) 



Figure 1: Direct and crossed two-photon diagrams of exchange interaction 

where W 2 = fi 2 R (Za) 2 /n 2 . The divergence of this integral identical to that for (10). Using 
Feynman's parameterization and the substitution of the type (10) to calculate (26), we 
obtain: 

Afi6 = jr^al n _ 2Jn2 _ 2_ _ 2_ + _ \ 

12m 2 \n 6 n 6 rr n 6 J 

3 Two-photon exchange interaction 

The amplitude of two-photon exchange interaction is represented on two diagrams of 
Fig.l. The appropriate quasipotentials are: 

V^Of ft = ^ Za)2 f /i(fc,mi,m 2 )rf 4 fc 

2 -> 

fi(k, mi,m 2 ) = m 2 (4m l + 2k ) - 2m x k< ) - 2k\ + -k 2 , 
D^-k) = k 2 - 2E 2 k + b 2 + ie ^ -2m 2 k + ie. 

V (b) (f) ft = ^ Za ^ f / 2 (fc,mi,m 2 )(f 4 fc 

27 IP, 9; n2 J ^ k _ p y + ie ^ k _ q y + ie } De ( k } D ^ p _ q _ k y y ) 

f 2 (k,mi,m 2 ) = m 2 (4mi + 2k°) - 2mik - 6/cq + ^-pk + ^-qk + ~^-k 2 , 

o o o 

D^p -q-k) = k 2 + 2E 2 k + 2k(p + q) - (p + q) 2 + b 2 + ie « 2m 2 k + ie. 

The main contribution of V 2l in the energy spectrum is proportional to a 5 . The 
analysis of corrections O(Za) 6 show, that they may exist in the energy levels also, if we 
take into account, for example, the contribution of photon poles, when k ~ a, \p\ ~ a, 
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Figure 2: Feynman diagrams of three-photon exchange interaction in muonium 



|g| ~ a, \k\ ~ a. In order to extract such terms from V27 let transform the product 
of electron and muon denominators in the direct two-photon diagram in the following 
manner: 



1 



-2ni6(k ) 1 
D e (k)D^-k) ~ -2E(P -b 2 )~2E 



1 



+ 



1 



(ko + ie)D e (k) (-ko + iejD^-k) 



(30) 



where the addendum with 5(k ) in the right part of (30) cancels with the iteration term 
of quasipotential. First addendum in the square brackets has the same structure in 
the leading order over l/rri2 as a crossed amplitude. Then the two-photon interaction 
quasipotential, which generates the correction (Za) 6 mf/m2 is equal: 



2i(Zaf 
3tt 2 



d 4 k[Ak 2 + 5k{p + q) - Qk 2 } 



\(k - p) 2 + ie\ [(k - q) 2 + te}D e (k)(2m 2 k + it) ' 



(31) 



We have calculated the contribution of V<j 7 to the energy spectrum for n=l and n=2 
by means of system for analytical calculations "Mathematica" [Q (package feynpar.m). 
The corresponding results are: 



AB 7 



-S(^) 

31m 
16m2 



n 



11 



1 



(32) 



4 Three-photon exchange interaction 

There are six diagrams, shown on Fig.2, which determine the three-photon exchange 
interaction in the muonium. 

Let consider the first diagram of Fig.2. The corresponding amplitude already has the 
factor a 6 , which appears due to electromagnetic vertices and Coulomb wave function. So, 
in the first stage of our calculations we have neglected by the electron and muon vector 
momentum of relative motion in the initial and final states, taking into account that the 



8 



necessary accuracy is already achieved. Then the first diagram amplitude (Fig. 2) takes 
the form: 



(Zaf 



I 




< 7i (?i ~p'+ mhiipi -P + mi)^ > 



(33) 



4tt 5 



(p 2 — w 2 + ie){p' 2 — w 2 + ie)[(p — p') 2 + ie] 



< 72 (P2 + g + ^"2)72 (g2 + g + ^"2)72 > 
D e (p)L> e (p')^(-p)^(-pO 



where D e>ll (p) are the denominators of electron and muon propagators: 

_D(±p) = p 2 — -u; 2 ± 2mp° + ie, w 2 = — b 2 , 



(34) 



and the corner brackets designate the averaging over Dirac bispinors; pi,P2 are the four- 
momenta of particles in the initial state; q±, q-i are the particle four-momentum in the 
final state. As usually is, the factor Za emphasizes exchanging character of photon 
interaction between particles. The exchange photon propagators were taken in Feynman 
covariant gauge. As it famous, using of Coulomb gauge is the most natural for exchange 
photons, because the Coulomb interaction dominates in the system (e~fi + ). Nevertheless, 
the equivalence of Coulomb and Feynman gauges in the scattering approximation for the 
three-photon diagrams calculations was shown in f|Tj]]. To construct the quasipotential of 
the system (e~/i + ) with L=0 and J=l, that corresponds to T-f 7 , let introduce the projector 
operator for initial and final states of the kind (6), setting also p = q = 0. Projecting the 
particles on 3 Si- state by means of (6), we avoid cumbersome matrix multiplication in 
the bispinor averages and immediately pass on to calculation of total trace in (33). As a 
result, the quasipotential of first diagram may be written in the form: 



D J (p) = p 2 — w 2 + ie, 

Fi(p,p') = fi2(p,p')m 2 2 + o/n m 2> fi2 = PP - 4m 2 - 1m x p^ - 2mip' Q - 2p p' , (36) 



We kept in (35) only the terms proportional to m\ and 7712, taking in mind the deter- 
mination of contribution to muonium fine structure in the leading order over parameter 
777.1/777.2. As will soon become evident, we can't restrict in F(p,p') only by terms ~ m 2 ,. 
The quasipotentials of the rest amplitudes of Fig. 2 may be constructed in a similar way. 
They differ from each other due to momentum dependence in muonic denominators and 
to the kind of functions fn (i=l,...,6). 

The parts of Fi(p,p'), proportional to 777.3, coincide in all six amplitudes. Let remark, 
that when substitute e — > 75 in projector operator (6) ( 1 S'o state), we obtain the same 
function fi2(p,p'), as for the 3 Si muonium. This means, that the muonium hyperfine 
splitting appears as effect of higher order on viiijm^- Functions fn are equal to: 



f /2 /2 / / / / 2/2 2/\ 

f 21 = — lOmip — 5poP + lOmipp + Ap pp — ^p pp + 1m\p + Ap p + 12p ^i— (37) 




D^D^D^p - pi)D e (p)D e (pi)D^-p)D^-p>) 



(35) 



fn(p,p') = 2mip /2 + p Q p' 2 + 10m ipp' + 2p pp' + 2p' pp'+ 

2 / 2 2 2 / 2 /2 / 

2mip + p Q p + 6777.-^0 + dttt^Pq + ^mip Q + 4mip — 4m 1 77 77 . 
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-6m 2 p(, + Amipl + 8mip Po - Smtffi + Ap p 2 , 

fsi = 2m\p' 2 + 4:p p' 2 + lOmipp' — Ap pp' + Ap'^pp' — lOrriip 2 — 5p' p 2 — (38) 

-Qmjpo + I2m\p' - Sm^l + 8mip p + + 4p Po 2 , 

j4i = 2mip + PoP —2m\pp + Ap pp + 2p pp — lOmip — 8p p — 12m 1 po+ (39) 

+6m 2 p() + 4mipg - 4mip Po + 4toiPo _ SPoPo; 

/si = -10m lP ' 2 - 8p p' 2 - 2m lH »' + 2p pp' + 4^pp' + 2m lP 2 + p(,p 2 + 6m 2 p - (40) 

-12m 2 p(, + Amipl - 4mip Po + 4raiPo _ 8 PoPo, 

/ 6 i = -lOmip' 2 - 5p p' 2 - 2mipp' - Ap pp - 4p(,pp' - 10m lP 2 - 5p' p 2 - (41) 

-6m 2 p - 6m 2 p[, - 8mipQ - 4mip Po ~ ^ m iPo- 

Integral function in (35) has simple poles over loop energies po,p'o in electron, muon 
and photon propagators. So, the most natural way of integration (35) consists in the 
calculation of integrals on po,p' at the first step, using the method of residues. But such 
an approach of calculation leads, nevertheless, to rather complicated intermediate expres- 
sions, what makes highly questionable its subsequent analytical integration on spatial 
momenta p, p*. So, we have used different approach of integration in (4), connected with 
transformation of muonic denominators, accounting the necessary calculational accuracy 
on parameter m\jm 2 . Considering that the spatial momentum of muonic motion in the 
intermediate state \p\ < m 2 , we obtain: 

2m 



D^ip) = p 2 -w 2 + 2m 2 po w 2m 2 p - \-ie) ^ 2m 2 (p + ie), (42) 



where the second approximate equality means that we neglect by the muon kinetic energy 
in the intermediate state. Doing so, we suppose that the integration contour on variable 
Po must be closed in the lower halfplane. Considering the terms, proportional to m\ in 
the numerators of all six diagrams (function fn(p,p')), we have arrived to the need of 
expression transformation, which includes the sum of muonic denominators (34). Using 
the second approximate equality from (42), we obtain: 

1 1 1 

+ — 7 rr^ — 7 — ; r + — 7 ttt^ w+ 



D^D^-p') D„(-p)D^p> - p) D^-p')D,{p - p>) 
1 1 1 



D^p)D^p - p>) D^pi)D^p> - p) Dp{pr)Dp(p) 
_ (-27ri)*(pd) (-2ni)5(p> ) 



(43) 



2m 2 2m 2 

In the energy spectrum the expression (43) will cause the corrections of order 0(a 4 ), 
which are canceled by the similar terms from the iteration part of the quasipotential. Con- 
sequently, to find the necessary contribution of order of a 6 , we must use first approximate 
equality in (42). Taking the difference 

1 1 (f + w 2 ) 



2m 2 (po - + ie) 2m 2 (p + ie) 4m 2 (p + ie) 2 ' 



(44) 



10 



let represent the quantity l/D At (p) in the form: 



+ 



(p 2 + w 2 ) 



(p) 2m 2 (p + *e) 4m 2 . (p + ie) 2 ' 



(45) 



Second addendum of (45) is of higher order on mi/m 2 in comparison with the first. But 
it leads to the necessary order correction on the other parameter a. Using the splitting 
(45), in the sum (43), let extract the terms, which generate the correction 0(a 6 ) and 
0(a G In a) in the energy spectrum. We may write them in the following manner: 



(ft 2 + w 2 ) 
8m% 

(p 2 + w 2 ) 
8ml 

(p-ft) 2 



2mS{p' ) 2mS(p' -p ) 2ixi5(p ) 



(p + ie) 2 (p + ief 



(p' + it 



+ 



(46) 



-2m8{p' Q ) 27tz5(pq — p ) 2iri5(po 



+ w 



(p + ie) 2 
2 m 5(p' ) 



•Sml 



+ 



(p + ie) 2 
2ni5(p' Q 



+ 



Po 



+ 



(Po + ie ) 2 
) 2%i5(po) 



(po + ie) 2 (po + ie) 2 (p' Q + ie) 2 _ 

It is evident from three-photon interaction amplitude of the type (33), that the parts 
of (46) give the necessary order corrections on a in the studied fine structure intervals. 
The same order corrections 0(mi/m2), as well as (18), will arise from the quasipotential 
terms containing the functions fn{p,p'), when we use the second approximation of (42) 
for muonic denominators. To do definite conclusion about the order of appearing terms 
in energy spectrum, which are determined by these quasipotential addenda, let transform 
them for greater simplification. Let consider for definiteness massless terms in the function 
fn(p,p f ), proportional to ~p 2 y p' 2 y pp': 



— 6pp' 



3p 2 

+3p 
1 



+ 



1 



1 



1 



D,{p) D M (-p) D^p-p') Dptf-p) 



+ 



(47) 
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+ 



D,(-P') D^{p') DM-P) D^p-p') 



+ 



dMp') dm) 

3p 2 



+ 



[-2^(po) + 27rz5(po-p[ ) )] + 

Lm2 ^772-2 



dMp) dm) dm-p) dm-p*) 

27n5(p' ) + 2iri8(p - p' )] • 



3p' 2 



6pp' 
2m? 



[-2m8(p ) - 2m5{p' ) + 2m5{p - p' )} . 



The transformation of other terms in functions fn(p,p') may be carried out by analogy. 
The next period of calculation consists in the integration over four-momenta in expres- 
sions (46)- (47). The typical two- loop integral, that results on this way has the following 
structure (ITi 



Ki = (4vr) 



d^pd A p' 



G i (p' ,Po,m 1 )P(p,p',w) 



-(2tt) 8 (p /2 - w 2 + ie)[(p -p') 2 + ie](p 2 - w 2 + ie)D e (j/)D e (p) 



(4f 
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where Gi(p' , po, mi) contains one 5 - function, and P(jf, p, w) is a polynomial. To calculate 
fundamental integrals (48) we have used Feynman parameterization in order to combine 
the denominators of the particle propagators, and the symmetry properties of the integral 
with the replacement p p'. There is the next set of functions Gi(p' ,po,mi), which 
appear in this paper: 

2m8(p -p' )2m 1 27ii5{p )2 mi 
Gi = (po + ieT » G 2 = - {p , +te)2 , G 3 = -2m5(p )2 mi , (49) 

G 4 = —27ii5(p' )2m 1 , G 5 = —2ni5(p — p' Q )2mi. 
The results of the integrations for K (48) are presented in the table [O. 



Table of the integrations i<Q (48), appearing in the 
muonium fine structure calculations 





p^(pff) 




ff 2 (pff) 


w 2 {pff) 




2 In 2 - i 


21n2 - i 


21n2 - | 







p*(f? 2 — pff) 


pp 1 (ff 2 — pp*) 


w 2 {jf 2 — f)ff) 


2 -6 

w z p z 


K 2 


1 l mi I 

2 2m 32 


1 | mi 13 

4 2w 32 


5 

32 


2 

3 




^2 


pff 






K 3 




I I mi 1 

4 2™ 4 


1 1 mi I 

2 2w 8 


1 


K A 


1 ]„ mi I 

2 2«) 8 


1 |„ mi 1 

4 4 




8 


K 5 


In 2 


ln2-i 


In 2 






Then the contributions, defined by expressions (37-41) and (46) are correspondingly 
equal: 

AB, = -^W" (50) 

2 777-2 
777^ 1 1 

AB 2 = (Za)^(Q\n2--) (51) 
m2 48 

We have introduced in (48) the photon mass w to avoid "infrared" singularities. The 
"infrared" logarithms Inw, containing this photon mass (see table of integrals Ki), and 
appearing at intermediate expressions, are mutually cancelled in the corrections ABi, 
AB 2 . 

Let consider now the quasipotential addenda, containing the momenta of particle rel- 
ative motion in the initial and final states. We denote them by r\ and r 2 correspondingly. 
Their consideration leads to modification of fn, which acquire the following additional 
terms: 

A/21 = 10m 1 pV 2 + 5pop'r 2 + m x pr 2 + 3p' pr 2 , (52) 
A/31 = —mip'ri — 3p p' r i ~ lOmipri — 5p' pri, (53) 
A/41 = mi(7pr 1 + 5p'r 2 + lOpri + 5pr 2 ) + Gpop'n + 5p pr 2 + 8p' pr 1 + 5p' pr 2 , (54) 
A/51 = rn\(— hp'r\ — 10p'r 2 — 5pr x — 7pr 2 ) — bpop'ri _ %PoP' r 2 — 5p' Q pr 1 — 6p' pr 2 , (55) 
A/6i = mi(— llpVi — \\pr 2 — \\pt\ — llpr 2 ) — 8pop' r i — &PoP r 2 — Sp'oPri — Sp' pr 2 . (56) 
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Using again the symmetry properties of appearing integrals under simultaneous variable 
replacement p p', r 1 r 2 , we obtain cancellation of all integrations in (52)-(56). So, 
the contribution of the particle relative motion in the fine structure with the accuracy 
0(777,1/777,2) is equal to zero. Thus the full value of the calculated correction (Za) 6 mf/m 2 
for hydrogen-like system S-states is defined as a sum of expressions (50) and (51): 

Afi 8 = (^) 6 l^(61n2-|). (57) 

77717772 V 48/ 

5 Discussion of the results 

In this paper we have calculated all possible corrections of order (Za) Q m\/m2 in the fine 
structure of hydrogen-like system on the basis of diagrammatic approach. Our total result 
is determined by the sum of terms AB { (4), (8), (11), (13), (24), (27), (32) and (57): 

AB tot = — In 2 H — — ^- + e n (58) 

tot U2 144 677 36n 2 96n 3 J nion 3 n y 1 



23mf(Za) 6 

12m 2 
31mj(Za) 6 



71=1, 



128m 2 ' 



n 



Numerical value of obtained contribution (58) for the "large" muonium fine structure 
interval 2 3 S' 1 -j- l 3 5*i is equal to 0,19 MHz. Recently, the calculation of the recoil correc- 
tions 0(Za) 6 m1/m 2 for the S-levels of hydrogen atom was done in || [7|, ||. The total 
contribution of the necessary order in the energy spectrum, which was obtained by using 
Braun formula, is equal ||: 

A£ « = fI + A_J_ + JL )(^M + f 41n2 _Z)(£^M. (59 ) 

\8 877 77^ 2n 6 J 7772 V 2/ 7772 

Comparing results (58) and (59), we see that the analytical expression of our correction 
(58) differs from (59), because we have calculated the contribution of necessary order to 
the energy spectrum, coming both from spin-independent and spin-dependant parts of the 
quasipotential. Numerical values of the contributions (58) and (59) for the levels with n=l 
and n=2 are equal correspondingly: n=l: -0,212 MHz, -0,065 MHz; n=2: -0,021 MHz, 
-0,006 MHz. The contribution of correction (58) to the fine structure interval 2S -j- 15* of 
hydrogen atom is equal to 21,5 KHz, whereas the values of similar contributions, obtained 
in [|| [J are equal correspondingly: 6,6 KHz ||; 14,5 KHz |J. 
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